Dynamic quantum Kerr effect in circuit quantum electrodynamics 
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A superconducting qubit coupled to a microwave resonator provides a controllable system that 
enables fundamental studies of light-matter interactions. In the dispersive regime, photons in the 
resonator exhibit induced frequency and phase shifts which are revealed in the resonator transmis- 
sion spectrum measured with fixed qubit-resonator detuning. In this static detuning scheme, the 
phase shift is measured in the far-detuned, linear dispersion regime to avoid measurement-induced 
demolition of the qubit quantum state. Here we explore the qubit-resonator dispersive interaction 
over a much broader range of detunings, by using a dynamic procedure where the qubit transition is 
driven adiabatically. We use resonator Wigner tomography to monitor the interaction, revealing ex- 
otic non- linear effects on different photon states, e.g., Fock states, coherent states, and Schrodinger 
cat states, thereby demonstrating a quantum Kerr effect in the dynamic framework. 

PACS numbers: 42.50.Pq, 42.50.-p, 85.25.-j, 03.65.Wj 



I. INTRODUCTION 

A major focus of quantum optics is the study of atom- 
photon interactions at the microscopic leveli Atomic 
and solid-state cavity quantum electrodynamics (QED) 
systems^ — with strongly enhanced coupling strength be- 
tween confined atoms and photons, have been developed 
to pursue this goal. The superconducting version, cir- 
cuit QED, provides a test bed for quantum microwave 
photons^ interacting with superconducting qubit s.^^ By 
placing a one-dimensional transmission line resonator in 
close proximity to a superconducting qubit ^ — the elec- 
tromagnetic fields in the resonator interact strongly with 
the qubit, forming a composite system with a modified 
energy spectrum (see Fig. 1(b)). The resonator pho- 
ton energy is effectively dispersed by the qubit, with a 
dispersion strength that depends strongly on the qubit- 
resonator frequency detuning. In previous experiments, 
the frequency and phase shift of transmitted photons 
were measured with a fixed qubit-resonator detuning, 
which we term the static scheme. Although the photon 
dispersion is well understood in this regime and provides 
a direct qubit readout scheme^ it is restricted to mea- 
surements with the qubit far detuned from the resonator 
to avoid measurement-induced demolition of the qubit 
quantum state. ^- 

In this paper, we introduce an alternative dynamic 
scheme for measuring qubit-photon interactions. As we 
adiabatically change the qubit frequency, we measure the 
accumulated phase shift of the photons while minimizing 
energetic interactions with the qubit. The detuning can 
range from large detuning corresponding to the linear re- 
sponse regime, to strong nonlinear dispersion found near 
and at resonance. The photonic response is measured 
with Wigner tomography, which enables a complete map 
of the photon state after a dispersive interaction with 



the qubit We further investigate the response of com- 
pletely nonclassical photon states, beyond the scope of 
prior static measurements, which have been restricted to 
measurements of classical coherent states. Measurements 
of the nonlinear dispersive response of superpositions of 
Fock states and Schrodinger cat states exhibit an excel- 
lent correspondence with expected dynamic response and 
demonstrate a quantum version of the Kerr effect. Clas- 
sically, the Kerr effect, a nonlinear response of macro- 
scopic materials is where the refractive index of a medium 
changes when varying the electromagnetic field strength. 
In this quantum version, a single qubit medium induces 
Fock-number dependent phase shifts to quantum state 
photons. Although the near-resonant nonlinear response 
of coupled atom-photon systems was previously explored 
using two-tone^SiSi and high power spectroscopy^^*22ri24 
our approach delineates the practical aspects associated 
with the strong non-linear phase shift of photons in the 
time domain. In the future, it should be possible to ex- 
ploit the large photon phase shift resulting from small 
qubit-resonator detunings to improve the quantum non- 
demolition qubit readout with higher fidelity. 



II. EXPERIMENT AND RESULTS 

We use a half-wavelength superconducting coplanar 
waveguide (CPW) resonator coupled to a superconduct- 
ing Josephson phase qubit (Fig. 1(a)). The resonator 
has a fixed frequency uOr/2Ti ~ 6.32 GHz and an energy 
relaxation time Ti 2 fis. The flux-biased phase qubit 
is a strongly nonlinear electrical quantum circuit, with 
a tunable \g) ^ \e) transition angular frequency ujq and 
|e) ^ I/) transition frequency ^ 200 MHz smaller (|^), 
|e) and |/) are the ground, excited, and third level qubit 
states) The qubit angular frequency ujq can be tuned 
rapidly using quasi-DC pulses on the inductively- 
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coupled qubit flux bias line; qubit energy level transi- 
tions are driven by 7r-pulse {(Jx,y) resonant microwave 
signals applied to the same line i^^i^^ The qubit state can 
be measured destructively using a single-shot readout in- 
volving short flux pulses followed by measurement using 
an on-chip superconducting quantum interference device 
(SQUID). 2^ The qubit has a relaxation time Ti 450 ns 
and a phase decoherence time T2 ~ 150 ns. The qubit 
and resonator are capacitively coupled with a coupling 
strength g/27T ~ 9 MHz (i.e., the qubit-resonator gap 
from spectroscopy is ~ 18 MHz at resonance). 
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FIG. 1. (Color online) (a) Circuit diagram for the cou- 
pled qubit-resonator system. The superconducting phase 
qubit (red) has a ground to excited state transition frequency 
0Oq/27T ^ 6 GHz. The flux bias line (black) is coupled to the 
qubit through a mutual inductance. A three-junction super- 
conducting quantum interference device (SQUID) (grey) is 
optimized to read out the qubit state. The CPW resonator 
(blue) has a constant resonant frequency cJr/Svr ^ 6.32 GHz. 
The resonator is weakly coupled to both the qubit and exter- 
nal microwave source by small capacitors, (b) Energy levels 
of the coupled system. The solid (dashed) lines denote the 
energy levels with (without) the qubit-resonator interaction, 
with \n,q) representing the photon number |n) and qubit state 
1^). The red and blue arrows represent one- and two-photon 
transitions with the qubit ground state. The green arrows 
represent one-photon transitions with the qubit excited state, 
(c) Plot of the \l,g) and |0, e) energies as a function of detun- 
ing A = ojq—ojr, calculated from the Jaynes-Cummings model 
using our experimental parameters. In an adiabatic process, 
the system initially at \l,g) (red dot) follows the path indi- 
cated by the red arrow. The resonant avoided-level crossing is 
g/n ^ 18 MHz. (d) Schematic of the experimental sequence, 
describing system initialization (labeled by "I") , adiabatically 
tuning and detuning the qubit, and finally performing Wigner 
tomography measurement (labeled by "M") of the final photon 
state. A reset puls©^ to ground is applied to the qubit if it 
was initially prepared in the excited state, to allow the qubit 
to perform the measurement. The dispersive part of the pulse 
is also shown schematically in (c). 



As shown for the specific case of |0,e) ^ |1,5') in Fig- 
1(c) [where states are denoted by |n, q) for resonator state 



\n) and qubit state \q) (Fig. 1(b))], the dynamic disper- 
sion scheme involves an adiabatic procedure^^ to vary 
the detuning A = uoq — uOr and thus the eigenenergies 
of the qubit-resonator system. We initialize the coupled 
system in the far-detuned regime, where the dispersive 
interaction is extremely weak. We then smoothly adjust 
the qubit frequency as a function of time t, so that the 
coupled system adiabatically follows its instantaneous 
eigenstates with eigenenergies E(t). The energy of a sin- 
gle photon (relative to the vacuum) is thus dispersed to 
Cpit) = E\i^g^{t) — E\Q^g^{t)^ and the dynamic phase^^ ac- 
cumulates as 6>(r) = — /o^[ep(t) — ep(0)](it. Although the 
system is returned to its initial settings at the end of the 
interaction, an accumulated phase shift is acquired even 
though the net energy exchange is zero. A Berry's phase 
component is ruled out because the parameter space is 
one-dimensional. We note that the qubit-resonator cou- 
pling strength is to all practical purposes constant in our 
system and the detuning is time-dependent, in contrast 
to the time-varying coupling and constant detuning dis- 
played by Rydberg atoms passing through a microwave 
cavity.^^^^ 

We first investigated the phase shift of a one-photon 
Fock state due to its interaction with a qubit in its ground 
state. The resonator was initialized in the Fock state 
superposition |0) + |1), where |0) serves as a reference 
stated The detuning was set to A^/27r = -250 MHz, 
giving a small photonic component ^/A^ 4%. We 
then applied a trapezoidal tuning pulse to the qubit, with 
linear ramp rate k = b MHz/ns and interaction detun- 
ing Af for an adjustable time r (Fig. 1(d)). At the 
end of the tuning pulse, the qubit was effectively decou- 
pled from the resonator and in its ground state, with a 
small error due to the non-ideal adiabatic process. The 
Wigner tomogram of the final resonator state was recon- 
structed by displacing the resonator state with a coherent 
microwave pulse (with complex amplitude a) and per- 
forming a qubit-resonator swap, from which we extract 
the photon state probability P^(a)J^i^ With sufficient 
sampling points a in the resonator phase space, we cal- 
culate the resonator density matrix and reconstruct the 
quasi-probability W{a)]^ 

We plot the resulting Wigner functions for detuning 
Aj/27r = —57 MHz and five different durations r in 
the top row of Fig. 2(a). Negative quasi-probabilities 
are a signature of the expected non-classicality of the 
superposed Fock states. The diagrams are not rotation- 
ally symmetric, encoding the relative phase shift between 
the |0) and |1) states. As the dispersion duration r in- 
creases, the Wigner function rotates counter-clockwise 
about the origin. The rotation is due to an overall phase 
shift, induced by the qubit state during the dispersive 
interaction. We plot resonator density matrices in the 
Fock basis in the bottom row of Fig. 2(a). In each den- 
sity matrix, the population probabilities poo and pn are 
nearly identical, but energy relaxation becomes more pro- 
nounced as r increases. For r = ns, the off-diagonal 
terms pio and poi are represented by horizontal arrows. 
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indicating no phase shift as expected. As r increases, 
the arrow representing piQ rotates with the same angle 
as that of the overah Wigner function, exhibiting an ac- 
cumulated non-zero phase shift; the time dependence of 
this phase rotation angle is shown in Fig. 2(b). The 
constant rotation rate represents a phase rotation speed 
ro = dO/dr = S.Itt rad//is and is equal to the expected 
frequency shift of the one-photon state |1) at the detun- 
ing Af. 
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FIG. 2. (Color online) Wigner tomography and phase rota- 
tion measurement for the initial resonator state |0) + |n). (a) 
Reconstructed Wigner functions for |0) + |1) after a dispersive 
interaction with a ground state qubit at A//27r = —57 MHz 
for five different value of r . The corresponding density ma- 
trices are shown in the bottom row. Sampling points^ for 
VK(a) comprised 60 points around two concentric circles with 
radii \a\ = 1.10 and 1.45. (b) Rotation angle of pno as a func- 
tion of T, measured from its value at r = 100 ns. (Inset) The 
n-dependent phase rotation rate is obtained from a linear fit 
to the data in (b). The results from both the linear approxi- 
mation (black line) and the exact diagonalization (red line) of 
the Jaynes-Cummings Hamiltonian are displayed without any 
fit parameters, (c) Phase rotation speed of PnO as a function 
of both Af and n. Solid lines are theoretical predictions from 
the exact diagonalization of the Jaynes-Cummings Hamilto- 
nian. 



We next investigated higher number photon states by 
initializing the resonator in |0) + |n) with n from 2 to 
5. Following the same pulse sequence, the density matri- 
ces and the Wigner functions of the final photon states 
were measured. Figure 2(b) shows the phase rotation 
angles of the pno element of the density matrix versus 
r. Linear fits allow us to extract the n-dependent phase 
rotation speed. In the far-detuned regime, the angular 
frequency shift scales as expected with photon number 
as ng'^/\Af\ (cf. Eq. 3 in Appendix B). This linear ap- 
proximation (black line. Fig. 2(b) inset) deviates from 



experimental data as n increases, because the dressed 
photon energies become anharmonic; a rigorous solution 
of the Jaynes-Cummings model (cf. Eq. 1 in Appendix 
B) is in good agreement with the data (red line. Fig. 2(b) 
inset). The phase rotation rate as a function of detun- 
ing Af (Fig. 2(c)) is also in excellent agreement with the 
full Jaynes-Cummings model, extending the nonlinear re- 
sponse measurements to zero detuning. This rate corre- 
sponds to the refractive index of macroscopic materials, 
where a single qubit, as an effective electric dipole, is the 
index medium which interacts with the electric field of 
quantum state photons. The Fock state dependent phase 
shifts imply a photon number dependent change in the 
refractive index of the medium, demonstrating a quan- 
tum Kerr effect in the dynamic framework. Although 
the energy levels of the qubit-resonator system can be 
measured using spectroscopy,^^ here the nonlinear phase 
shift in the strong dispersive regime is directly and quan- 
titatively measured in the time domain in a circuit QED 
system. 
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FIG. 3. (Color online) (a) Wigner functions and density 
matrices for a coherent state (a ^ 2.0) interacting with a 
ground state qubit at three different durations r for detuning 
A//27r = —50 MHz. Axes at the origin are indicated. The 
center of the coherent-state Wigner function rotates counter- 
clockwise about the origin with time, (b) Wigner functions 
for a Schrodinger cat state \a = V^i) -\- \a = —\f2i) after 
different interaction times with a ground state qubit. A fixed 
detuning A/ from —56.3 MHz to —14.1 MHz was used for 
each column. As the duration r increases from 150 ns to 
250 ns, the Wigner function rotates counter-clockwise, with 
rotation rate increasing as the absolute value of detuning Af 
decreases. For small detunings, enhanced distortions appear 
in both the coherent peaks and the interference fringes. 



In a standard resonator transmission measurement, a 
coherent state \ol) is injected in the resonator by a clas- 
sical microwave pulse. Here we prepared a coherent res- 
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onator photon state with complex amphtude a ^ 2.0 
and apphed the dynamic scheme to measure the phase 
shift induced by interaction with a ground-state qubit. 
Figure 3(a) presents the measurement after a trapezoid 
pulse at A//27r = —50 MHz. The Wigner function of 
the initial state (r = ns) shows a coherent peak cen- 
tered at a and a symmetric noise of quantum uncer- 
tainty. As r increases, the coherent peak rotates counter- 
clockwise about the origin, with only a slight and gradual 
increase in distortion. For each Fock constituent, the lin- 
ear approximation gives a phase shift rate ~ nro = 
ng'^/\Af\. The coherent state is thus uniformly trans- 
formed to \ae'^^^^) after a time r, based on the Poisso- 
nian distribution. In the Wigner representation, each 
Fock constituent has a rotation rate rw,n = Tn/n = ro 
and the Wigner function rotates as a whole with rotation 
angle ror, which is equal to the phase shift measured in a 
resonator transmission measurement. This simple behav- 
ior is most accurate in the far-detuned regime, where the 
dressed resonator levels are still nearly equally spaced. 
At small detuning, where anharmonicity becomes more 
significant, the various Fock components begin to rotate 
out of phase and the Wigner function begins to display 
squeezing (see Fig. 3(a)). 

To further illustrate the linear-nonlinear crossover 
of the dispersive interaction, we investigated the 
Schrodinger cat state \a = V^i) -\- \a = —\/2i). In Fig. 
3(b) we display the reconstructed Wigner functions for 
four different values of detuning Aj, and three interac- 
tion times r. At zero interaction time, the Wigner func- 
tion contains two well-separated Gaussian peaks, arising 
from the two coherent state elements, with clear interfer- 
ence fringes between the two peaks. The Wigner function 
rotates counterclockwise with increasing time r, as be- 
fore. For detuning A//27r = —56 MHz, the shape of the 
Wigner function is almost preserved as r increases, con- 
sistent with expectations from the linear approximation. 
As the absolute value of detuning Aj decreases, the rota- 
tion speed increases and distortion becomes significant. 
The two Gaussian peaks are twisted and the interference 
fringes deform with time, due to the differing rotation 
speeds of the Fock components. This strongly nonlin- 
ear effect is most prominent for zero detuning, where the 
rotation rate scales as Vyj^n oc Xj ^Jnw^^ 

The examples above were measured with a qubit ini- 
tialized in its ground state, for which qubit relaxation 
and phase decoherence are less important. We also ex- 
amined the response of a resonator interacting with a 
qubit initialized in the excited state. Figure 4(a) shows 
the Wigner function for a resonator initialized in |0) + |1) 
after interacting with an excited-state qubit at detun- 
ing Aj/27r = —35 MHz. In contrast to the response 
seen with a ground- state qubit, we observe the opposite 
rotation direction, indicating a negative shift in photon 
frequency and phase. The nonclassical negative quasi- 
probabilities of the Wigner function diminish much more 
quickly, in accordance with the rapid decoherence of the 
off-diagonal terms in the density matrix. In addition, for 



three different final detunings A/ (Fig. 4(b)), the phase 
rotation angle of poi is no longer linearly dependent on 
r but shows a periodic ^-shaped structure (Fig. 4(b)), 
which can be related to qubit relaxation. 



(a)r- 



70 ns 



100 ns 



130 ns 160 ns W(a) 




(b) 



0123 0123 0123 0123 0123 

n 



■a 


O.Oti 




res: |0)+|1)- 


CO 






[jVPubit: |e) 


CX) 








CD 


-0.471 






CO 

cc 




45 MHz 




Q. 


-O.Stt 


. • 35MHz 








• 25MHz 














100 200 300 

Interaction time x (ns) 



-2-10 1 2 
Re(a) 



< -0.47T O.Ot: 0.47T 

Angular coordinate (rad) 



FIG. 4. (Color online) Measurement of photon states after 
dispersive interaction with an excited-state qubit. To mini- 
mize qubit relaxation during the interaction, we use a higher 
tuning ramp speed of A; = 7 MHz/ ns. The qubit was reset 
to its ground state prior to the tomography measurement, by 
transferring the excitation to a nearby two-level state using an 
z-swap gate.i^ (a) Wi gner functions and density matrices for 
the initial resonator state |0) + The dispersive interaction 
was measured for detuning A//27r = —35 MHz for five dif- 
ferent durations r. (b) Rotation angle for pio vs. r for three 
different detunings A/. Solid lines are numerical simulations 
of the Lindblad master equation including both system deco- 
herence and the third level of the qubit. However, we have 
to use a shorter qubit decoherence time to yield reasonable 
agreement between the experimental and theoretical results. 
For |A/|/27r of 25 MHz, 35 MHz, and 45 MHz, the qubit relax- 
ation times Ti were 200 ns, 220 ns, and 250 ns, respectively. 

(c) Wigner function for a coherent state after interacting with 
an excited-state qubit for r = 220 ns at A//27r — —75 MHz. 

(d) Radially averaged W[cx) as a function of the angular co- 
ordinate, extracted from the Wigner function in (c). 



To understand the effects of qubit relaxation, we in- 
vestigated the response of a coherent photon state (Fig. 
4(c)). During the dispersive interaction, the qubit can 
relax at any moment, performing a stochastic quantum 
jump in its time trajectory. As a consequence, the co- 
herent peak in the Wigner function will initially rotate 
clockwise, but when the qubit relaxes, switch to a coun- 
terclockwise rotation for the remainder of the dispersive 
interaction. The amplitude of the jump probability is de- 
termined by the qubit relaxation time T\ and decays ex- 
ponentially as time increases. The reconstructed Wigner 
function thus comprises the main coherent peak (from the 
excited-state qubit) accompanied by a small peak with a 
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tail (due to rotation from the ground-state qubit). Nu- 
merical simulations of the Lindblad master equation are 
quantitatively consistent with the experimental results 
(not shown). 

For a general resonator state interacting with an ex- 
cited qubit, the qubit relaxation and the mixed qubit 
state lead to a decoherence effect on the final photon 
state. The decoherence effect is averaged over many in- 
dividual quantum jumps in the ensemble measurement 
required for collecting quantum statistics. For a res- 
onator initialized in |0) + |1), the excited qubit disperses 
the photons with a phase shift, i.e., rotates the hori- 
zontal off-diagonal vector poi to an angle 0{t). The 
off-diagonal terms from qubit quantum jumps are rep- 
resented by vectors with rotation angles distributed over 
the range [^(r), —0{r)] and amplitudes corresponding to 
the probability of the specific quantum jump. Averaging 
over all the vectors resulting from quantum jumps de- 
creases the phase rotation angle of poi when is small, 
i.e. still an acute angle. When 6 is larger than 7r/2, the 
averaged vectors from the qubit relaxation increase the 
phase rotation angle. This explains the S-shape structure 
in the time-dependent phase rotation angles (Fig. 4(b)). 
These data are in fact predicted (Fig. 4(b)) from numer- 
ical simulations using the Lindbald master equation (cf. 
Appendix B) incorporating the third qubit level, qubit 
decoherence, and TLSs. 



Appendix A: EXPERIMENTAL SETUP 

The superconducting quantum circuit used three 
coplanar resonators interconnected by two superconduct- 
ing phase qubits^— although in this experiment we 
used only one qubit and one resonator. The fabrica- 
tion process is a combination of photolithograhpy and 
plasma etching on a multi-layered structure, with a su- 
perconducting rhenium base film, Al/AlOx/Al Josephson 
junctions and amorphous silicon for the shunt capacitor 
and wiring crossover dielectrics^ An aluminum sample 
box containing the wire-bonded device was mounted on 
the mixing chamber of a dilution refrigerator operating 
at 25 mK. 

The quasi-DC tuning pulse was generated by cus- 
tom electronics including an FPGA-controlled digital-to- 
analog (DAG) converter. Two similar DAG outputs, the 
I and Q signals of an IQ mixer, provided the sideband 
mixing, phase, and pulse-shape control of the carrier sig- 
nal from a microwave generator to produce the ax,y pulse. 
The qubit and resonator frequencies were determined by 
spectroscopic measurements. The qubit-resonator cou- 
pling strength g and the resonator microwave pulse am- 
plitude a were calibrated using a photon population anal- 
ysis for a coherent stated More details on the microwave 
calibration procedure are described in the supplementary 
information of Ref. 19. 



III. CONCLUSIONS 

In conclusion, we have used an adiabatic dynamic con- 
trol of the qubit to study its dispersive interaction with 
a cavity resonator. In this dynamic scheme, the fre- 
quency shift of multi-photon Fock states is extended to 
the strongly non-linear dispersive regime, close to zero 
detuning. We can fully control and measure the accumu- 
lated phase shift of complex photon states, demonstrated 
by the excellent agreement between the experimental 
data and the theoretical prediction. In the strongly non- 
linear dispersive regime, a quantum Kerr effect is ob- 
served for the coherent and other non-classical photon 
states: We illustrate the nonlinear response of photons 
to the refractive index effect of our artificial atom using 
Wigner tomography. Furthermore, we reveal the distinct 
phase shift of photons induced by the excited state qubit, 
and interpret the peculiar photon decoherence resulting 
from qubit state relaxation. 



ACKNOWLEDGMENTS 

This work was supported by lARPA under ARC award 
W911NF-08-01-0336 and under ARC award W911NF-09- 
1-0375. M.M. acknowledges support from an Elings Post- 
doctoral Fellowship. Devices were made at the UG Santa 
Barbara Nanofabrication facility, a part of the NSF- 
funded National Nanotechnology Infrastructure Network. 



Appendix B: MODELING AND SIMULATION 

The qubit-resonator coupled system is described by the 
Jaynes-Gummings Hamiltonian,^ H = hwr{a^a^ 1/2) + 
hwqCFz/2 + hg{(j~^a + a+a"), where hwr is the single pho- 
ton energy of the relevant electromagnetic field mode, 
hujq the excitation energy of the qubit, and a+(a) and 
cr+((j~) the creation (destruction) operators for the res- 
onator photons and the qubit, respectively. With a cou- 
pling strength ^, the interaction term describes energy 
exchange between the qubit and the resonator. The pho- 
ton energy is effectively dispersed by virtual energy ex- 
change when an electromagnetic wave transmits through 
the resonator at a detuning A = ujq — uj^. In the Jaynes- 
Gummings model, the qubit-resonator interaction leads 
to two eigenenergies nhur ± f ^g^n + for each state 
pair \n^g) and |n — 1, e) (A < is assumed). As shown 
in Fig. 1(b), the energy of a n-photon state is dispersed 
to 



= nhwr + J(VV^ + A2 _ |A|), (Bl) 

when interacting with a ground state qubit. With the 
excited state qubit, the dispersed energy becomes 



6 



nep^e — E\n,e) " -^|0,e) 

= nhwr - ^{V^aHn + 1) + A2 - ^4^2 + A2). 

(B2) 

In the far-detuned regime, the hnear approximation is 
apphed and the n-photon energy is simphfied to be 

nCp^cTz — nhujr + azuhg'^ / A (B3) 

for the ground-state and excited qubit (represented by 

The dispersive qubit-resonator interaction can be 
strongly affected by the qubit relaxation and decoher- 
ence. Our numerical simulations used the Markovian 
Lindblad master equation for system density matrix p, 

^ = -^[H{t),p]+ Y: L-p-L^-h^-L-p-lp-L^-L, 



where H{t) is the instantaneous total Hamiltonian 
and two superoperators, Li = l/^/Tla~ and L2 = 
l/v^cr^cr~, describe qubit relaxation and decoherence. 
To simulate the phase qubit, the third qubit level can be 
added in the numerical master equation. We find that 
for the photon state phase shift induced by a ground 
state qubit, the numerical results of the master equa- 
tion including decoherence and the third level contribu- 
tion are very close to the analytic solution of the Jaynes- 
Cummings Hamiltonian without dissipation, with differ- 
ences smaller than our experimental resolution. For pho- 
tons interacting with an excited qubit, in contrast, nu- 
merical simulations with both decoherence and the third 
level have to be applied when comparing experimental 
data with theoretical results (see Fig. 4(b)). 
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